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ABSTRACT  AND  INTRODUCTION 

Problems  involving  computation  of  the  probability  of  damage,  or  the  expected 
fraction  of  the  target  damaged,  occur  frequently  in  military  operations  research. 
These  problems  have  been  worked  on  by  many  people  over  a  long  period  of  time. 

This  study  is  intended  to  summarize  some  cf  the  more  important  results  and  to 
indicate  sources  of  information  on  other  problems.  No  attempt  has  been  made  to 
include  an  extensive  bibliography  or  to  track  down  the  first  man  who  derived  any 
particular  result. 

Although  many  of  the  results  were  known  long  ago  -  some  of  them  during  the 
19th  century  or  earlier  -  it  was  not  until  World  War  n  that  systematic  efforts  were 
made  to  obtain  answers  for  a  wide  variety  of  situations.  Recently  several  reviews 
have  been  made  of  this  material.  One  of  these  is  reference  (a),  which  does  not 
include  some  of  the  results  given  here.  On  the  other  hand,  reference  (a)  gives 
results  on  the  inverse  problem,  that  of  determining  the  probable  location  of  the 
aiming  point  from  the  locations  of  the  burst  points  or  from  the  damage  produced. 

The  inverse  problem  is  of  less  interest  than  tie  direct  problem;  the  results  are 
of  little  usefulness  in  most  damage  problems. 

SINGLE-SHOT  PROBABILITIES 

In  this  section  we  will  consider  the  problem  of  computing  the  probability  of 
killing  the  target,  or  of  inflicting  a  stated  degree  of  damage,  by  a  single  weapon 
(projectile,  bomb.  etc. ).  Because  of  its  frequent  occurrence  the  two-dimensional 
case  will  be  used  throughout  for  the  illustration  of  methods  and  statements  of 
results.  In  most  cases  the  manner  nn  which  the  results  would  be  adapted  to  one 
dimension  or  to  more  than  two  dimensions  should  be  evident.  The  derivation  must 
be  examined,  however,  in  every  instance  in  which  only  the  results  are  given. 

We  will  consider  both  point  and  area  targets,  a  point  target  being  one  whose 
dimensions  are  small  compared  with  the  "damage  radius "  of  the  weapon  relative 
to  the  target.  We  will  also  consider  two  types  of  "conditional  damage  function", 
that  is,  the  function  that  describes  the  probability  of  damaging  the  target  if  the 
weapon  detonates  at  a  given  point  relative  to  the  target.  These  two  functions  are 
usually  referred  to  as  the  cookie-cutter  function  and  the  Gaussian  function.  In 
the  former  case  the  probability  of  damage  is  I  within  some  area  about  the  point 
of  detonation  and  is  zero  outside  this  area.  In  the  latter  the  probability  of  damage 
varies  gradually  from  a  value  of  1  at  the  point  of  detonation,  decreasing  toward 
zero  as  the  distance  from  the  point  of  detonation  increases. 

These  are  some  of  the  principal  results: 

★  SINGLE  POINT  TARGET 

Consider  the  target  at  the  origin  of  a  rectangular  coordinate  system  in  a 
suitable  plane.  This  plane  might  be  the  horizontal  piane  or  a  plane  normal  to  the 
trajectory  of  a  weapon  which  would  pass  through  the  target.  Let  the  conditional 
probability  of  damage  to  the  target,  if  the  weapon  is  fired  or  released  on  a  trajectory 
through  the  point  (x,  y),  be  p^(x,  y).  Let  the  joint  probability  density  function  of 

x  and  y  be  f(x,  y).  The  joint  density  function  f(x,  y)  describes  the  density  of  shots 
or  bomb-drops  in  the  target  plane.  It  is  non-negative  and  is  so  chosen  that  the 
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integral  of  it  over  any  area  is  the  probability  that  the  weapon  will  fall  in  this  area. 
Then  the  probability  of  damage  to  the  target  from  a  single  weapon  chosen  at  random 
is 


■si 


Pd&  y)  f(*>  y)  dx  dy 


*00 


Some  cases  are  given  below: 

•  Circular  Cookie-Cutter  Damage  Function: 


Pd(*>  y)  = 


1  if  x2  +  y2  £R2 
0  otherwise 


For  this  function  the  probability  of  damage  is 

P=  JJ  f(x,  y)dx  dy, 
circle 

where  the  integral  is  to  be  taken  over  the  circle  of  radius  R  with  center  at  the  origin. 

*  Probability  of  damaging  single  point  target  with  a  single  shot  when  the  con- 
tHtionirflamage  function  is  the  circular  cookie-cutter  function  and  when  the 
distribution  of  shots  is  circular  normal  and  centered  on  the  target? 

More  precisely,  we  assume  that  x  and  y  are  independent  variables  having  normal 
distributions  with  means  zero  and  common  standard  derivation,  a  .  Then  f(x,  y) 
has  the  form 

f(x,  y)  =  - i-*-  exp  [  -  <x2  +  y2)/2  a2 

2  k  a1  1 

By  elementary  integration  we  find  that 
P  =  1  -  exp(-R2/2  a  2) 


■  Probability  of  damaging  single  point  target  with  a  single  shot  when  the  con¬ 
ditional  damage  function  is  the  circularcoolcie-cutter function  and  when  the 
distribution  of  shots  is  elliptical  normal  and  centered  on  the  target: 


This  case  is  the  same  as  the  one  preceding,  except  that  the  two  standard 
deviations  are  unequal.  This  slight  change  makes  the  integration  difficult.  If  we 
change  to  polar  coordinates  (r,  0)  and  integrate  with  respect  to  0  we  obtain  the 
following  result: 
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P  =  — o —  ^ y  e*P  £~r2  cos20/2  -  r2  siu20/2  O ^  j  rdrd0 


2r  O 


a 

L_  f 

<r  <r  „  J 


e  ar  r  dr 


x  y  o 
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I  r2 


cos  ^ 


d*.  (4>  =20) 


=  V’z  [ 


e  ar  IQ(br2)  r  ^ 


-i _  f 

1  °  2  'l 


e  au  I  (ba)  du 
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Here  <J  is  the  smaller  of  a  ,  a  and  <ro  is  the  larger  of  these.  Also, 
i  x  y  z 

_  _  1  ,  1  ,  l,^_l/l  lx 

a_T  +  b"T  (1T2  -  *cr^> 

and  Iq(x)  is  the  modified  Bessel  function  of  zero  order.  If  we  replace  the  Bessel 

function  by  its  series  expansion  and  integrate  term  by  term,  we  get  the  following 
result: 

„  fJl /0  „  2.k+l  . 

00  u  (R  /2  O .  )  *  i  i,  2i  2  2  i 

p  =  — V*  1-\\K  _ i _  V'  /-l /A\l(?\  t£l\  n 


£  H/4)‘<pq‘>  ( 

i=0 


£  (-D 


For  rases  in  which  R  is  small  compared  with  this  series  converges  rapidly. 

The  integral  of  the  elliptical  normal  distribution  over  a  circle,  as  derived 
above,  is  equal  to  the  integral  of  the  circular  normal  distribution  with  unit  standard 
deviation  over  the  ellipse  that  has  semi-axes  R/dx  and  . 

The  probability  of  damage  in  the  case  of  unequal  standard  deviations  can  be 
approximated  by 

P  =  1-  exp  (-R2/2<r2) 

where  a  is  an  appropriate  function  of  a  and  o  .  Functions  frequently  used  are 

the  geometric  mean  a  ,  the  arithmetic  mean  a  ,  and  the  root- mean- square  o  . 

g  a  s 
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where  #  is  the  normal  error  integral: 


'00 


■  Probability  of  damaging  single  point  target  with  a  single  shot  when  the  con¬ 
ditional  damage  inaction  is  toe  circniar  conlae-caaer  function  and  wjjea  tbe 
distribution  of  shots  is  elliptical  normal  and  x  and  y  are  correlated: 


We  assume  that  the  means  are  zero  but  x  and  y  are  not  independent.  Let  the 
covariance  be  0^.  This  is  equal  to  the  expected  value  *  of  xy  and  is  equal  to 

p  0  0  where  p  is  tbe  correlation  coefficient.  The  probability  of  damage  ta 

xy  x  y  xy  2  7 

this  case  can  be  obtained  from  the  above  results  by  replacing  0^  and  0O~  by 

and  respectively,  where  and  are  the  roots  of  the  equation: 
(dx2-X)(0y2-X)-  0^=  s 


witls  X^  —  X,.  The  roots  of  this  equation  axe  often  called  die  eigen-values  of  the 

covariance  matrix.  The  equation  is  obtained  by  equating  to  zero  the  determinant 
in  which  the  elements  are  die  corresponding  elements  cf  the  covariance  matrix 
with  the  elements  of  die  main  diagonal  reduced  by  X. 


This  result  can  be  obtained  easily  from  die  joint  density  function,  which  is: 
*£x,y)  =  (1/2  t  0^  0y  Vl  -p2)  . 

expf-  ( cr^x  '  2p  0X  0y  xy  +  dx272)/2  0x2  0„2  (L  p2)  ] 


where 


0  “  p  0  0  . 

xy  x  y 


We  now  rotate  the  axes  through  one  of  the  angles  dig? 


eliminates  tbe  cross-product  term  in  the  quadratic  form.  The  transformation  is 
linear  with  determinant  equal  to  1 .  If  u  and  v  are  the  new  rectangular  coordinates, 
the  function  f(x,  y)  reduces  to 


(1/2  t  exp  (-  uZ/2  X.  -  v2/2 


where  X^  and  are  the  roots  of  tbe  equation  given  above.  Since  the  transformation 

dees  not  alter  the  conditional  damage  function,  the  results  can  be  obtained  from  the 

2  2 

previous  case  by  replacing  <?  and  &  "  by  X^  and  X^. 


*  Decause  the  means  are  zero.  More  generally,  tbe  covariance  Is  equal  to  tbe 
expected  value  of  (x  -  x)  (y  -  y),  where  x  and  y  are  the  means  of  x  andy  respectively. 
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r  =  a2/  l  f»2  *  2  0^}  (a2  *  2  n^2) 
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pom  target  with  a  siagle  shot  when  the  coo- 


otthe  drcMt  6wt«ia  haetka  and  ufaen  tbe  dia- 

j i  stow  5  cxrcaiar  wareaal  aad  is  centered  at  a  point  h  units  £r>m 


targes: 


.  t%e  atteyal  caa  be  computed  readily  by  tbe  separation  of  variables. 
Tbe  pesfcndsSaty  ss  gswer  by 


S'  =  j  a2/f*2  -  2  P2)  j  eap  -fc2/(a2  *  2  *2)  j 


Erarff  she  vtac  ftpteme  is  obci:aed  for  tbe  case  in  which  tbe  distribution  is  tbe 
Carton  dsaerabfiCKW  wstifc  dnsxr  fuactioe 


S*.  y)  --  a/2  X02)  exp  F  -(x2  ♦  y2  +  h2y2c2  J  I^b  V x2  +  y  1/o1  j  . 

•  ?toba5i£ary  o£  dassagigg  siagjje  point  target  with  a  single  shot  when  the  coa- 
.fcottal  daaaawe  fcacbow  is  the  circalar  Gaussian  function  and  when  the  dis- 
rafemnoo:  ST&ss  is  elliptical  normal  cot  centered  on  the  target: 

Asaoaae  sfeas  dhe  center  of  tbe  distribution  is  at  (b^,  h^).  Again  tbe  integral 
caa  be  sepaiaced.  sad  tbe  probability  of  damage  is 


?  =  |  a2/  k(a2  *  2  dw2)  (a2  r  2  <ry2)  j  • 


I  -h  ^(a*  +  2  d  2)  -b  2/(a2  +  2  o  2)  1 
exp  ^  x  x7  y  ,x  y  j 


•  Exiipocal  Gaussian  Damage  Function: 
?d 


(*.  T)  =  «p  (*  x2/ax2  -  y2/ay2) 


fcs  this  case  tbe  integral  of  tbe  conditional  damage  function  over  the  entire  plane 

is  v  a  a  ,  use  area  oi  an  ellipse  with  semi-axes  a  and  a  .  We  can  determine 
x  f  x  y 

tbe  probabilities  in  this  case  from  tiie  corresi>onding  probabilities  for  the  circular 

Gaussian  damage  function  by  making  the  following  substitutions: 
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Replace  a 
Replace  ia 
a2  +2 


a2  +2 


2 

a  +2 


in  the  numerator  by  a^  ay 
the  denominators: 

°x  b*  \  +  2  *x 

°y  by  Ay2  +  2  <ry2 

<r2  by  V(a2  +  2  a1)  (a  2  +  2  o2) 
x  y 


•  Other  Conditional  Damage  Functions: 

Por  some  weapons  and  weapons  effects  the  conditional  damage  function  appears 
to  lie  somewhere  between  the  circular  cookie-cutter  function  and  the  circular 
Gaussian  function.  In  appendix  A,  a  sequence  of  functions  of  this  type  is  proposed. 
The  first  member  of  this  sequence  is  the  Gaussian  function,  and  the  limiting 
member  is  the  cookie-cutter  function.  These  functions  have  the  property  that  the 
probability  integrals  for  the  usual  cases  can  be  computed  in  finite  form  in  terms 
of  known  functions.  The  second  member  of  the  sequence  yields  probabilities  that 
are  approximately  halfway  between  the  corresponding  probabilities  for  the  Gaussian 
and  cookie-cutter  functions.  This  property  could  be  used  as  a  means  of  approxi¬ 
mating  the  difficult  cases  for  the  cookie-cutter  function.  However,  since  die 
second  member  of  the  sequence  probably  is  closer  to  the  true  conditional  damage 
function  than  is  the  cookie-cutter  function,  it  is  more  reasonable  to  use  the  prob¬ 
abilities  given  by  this  second  member  of  the  sequence. 

★  SEVERAL  POINT  TARGETS 

When  several  point  targets  exist,  there  is  generally  a  need  to  answer  questions 
of  the  following  type:  What  is  the  probability  of  damaging  all  the  targets,  or  at 
least  a  specified  number  of  them,  with  a  single  weapon?  (The  case  of  several 
weapons  will  be  discussed  later.)  To  answer,  we  must  first  compute  the  desired 
probability  under  the  condition  that  the  weapon  detonates  at  a  particular  point 
(x,  y).  This  quantity  is  then  treated  as  the  conditional  probability  in  the  previous 
work  and  we  must  average  this  function  over  the  shot  distribution. 


For  the  cookie-cutter  damage  function  the  integrals  involved  are  invariably 
difficult.  The  integrals  are  fairly  easy  to  compute,  however,  in  the  case  of  the 
Gaussian  damage  function.  We  give  one  example  below. 


Assume  that  the  distribution  of  shots  is  circular  normal  of  standard  deviation 
a  and  is  centered  at  a  point  midway  between  two  targets.  Assume  that  the  targets 
are  located  at  the  points  (-h,  0)  and  (h,  0).  If  the  weapon  detonates  at  the  point 
(x,  y)  the  probability  of  damaging  both  targets  is 


Pd(x,  y) 


exp 

exp 


-[(x-h)2  +  y2  +(x  +  h)2 
-2  (x2  +  y2  +  h2)/ a2  ] 


+ 


1 
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To  obtain  the  desired  probability  we  must  now  multiply  by  the  density  function  and 
integrate  over  the  entire  plane.  The  result: 

P  (both)  =  ^  a2/ (a2  +  4 a2)  exp  (-2h2/a2) 

★  AREA  TARGET 

In  this  case  we  assume  that  the  target  is  spread  over  a  large  area,  the  maximum 
dimension  of  this  area  being  at  least  as  large  as  the  damage  radius  of  the  weapon. 

It  is  assumed  that,  for  a  large  fraction  of  the  detonations,  some  parts  of  the  target 
would  be  damaged  critically  while  other  parts  would  be  left  undamaged.  In  this 
case  it  makes  little  sense  to  talk  about  probability  of  damage.  Instead  we  must 
consider  the  problem  of  determining  the  expected  fraction  of  the  target  that  will 
be  damaged,  or  the  probability  of  damaging  at  least  a  given  fraction  of  the  target. 

Assume  that  the  burst  occurs  at  the  point  (x,  y).  Then  if  (X  Y)  is  any  point  of 
the  target,  the  conditional  probability  of  inflicting  damage  to  this  point  is 

Pd  (x  -  X,  y  -  Y) 

Then  the  conditional  fraction  p(x,  y)  damaged  if  the  burst  occurs  at  (x.  y)  is 
p(x,y)  =  (l/A)  Upd(x-X  y- Y)dXdY 


where  the  integral  is  taken  over  the  target  of  area  A.  The  expected  fraction  of  the 
target  damaged  is 


00 

F  =  j*  J  p<x,  y)  f  (X,  y)  dx  dy 
-  oo 

where,  as  before,  f(x,  y)  is  the  joint  density  function. 
Let  Sq  be  the  Sv  t  of  points 

Sq:  Paints  (x,  y)  such  that  p  (x,  y)  ^  Fq 


Then  the  probability  of  damaging  at  least  the  fraction  F  is  equal  to 

c 

p(2 F„>  =  j  j"  y)  ‘tody 

So 

where  the  integral  is  taken  over  the  set  SQ.  Sandia  Corporation  literature  describes 

an  analogue  computer  which  vas  designed  to  obtain  this  function  for  a  wide  variety  of 
target  complexes  and  density  functions. 
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We  illustrate  these  ideas  with  one  example.  Let  the  target  be  a  circle  of  radius 
T  and  let  the  conditional  damage  function  be  the  circular  Gaussian  function  of  equivalent 
radius  a.  Assume  that  the  distribution  is  circular  norir.2!  of  standard  deviation  a 
and  centered  at  the  ceoter  of  the  target.  Then  the  conditional  fraction  of  the  target 
damaged  is  J 

p(x.y)  =  f l/t  ^  J  3  exp  |  -  j  (r  cos  0  -  x)“  *  (r  sin  0  -  y)2  J  /a2  }  rtirtiG 

o  o  >/a2 

-  (a2/T^)  exp  -  <xx  +  y*J/a2  j  J  e  u  1(2  }  u(x2  +  y2J^a)  du 
*-  J  o 

The  expected  fraction  of  the  targe:  damaged  is 


r  - 


(1/2  r  a2)  \  \  p (x.  y)  axp  !  -  (x2  +  y\'2a2  1  dx  dy 

-J 


=  (»2/T2)  j  1  -  exp  f  -  T2 /(a2  +  2  <y2)  ]  j 


r  —  l  j  j 

To  obtain  the  probability  of  damaging  at  least  the  fraction  Fq  we  would  have  to 
solve  die  equation  obtained  by  equating  p(x.  y)  to  F  .  This  will  be  a  circle  whose 

D 

radius  could  be  obtained  by  trial  and  error.  Let  toe  radius  of  this  circle  he  R  . 

o 

Then  die  probability  of  aaisagmg  ex  least  this  fraction  *oald  be  equal  to 

2 

PtiFo)=  fl-«P<-  Ro2  nc\  [i-«p<-tV)  | 

[o  .  F0  >  (aV2)  [  I-exp  <-tV)  j  . 

If  we  replace  the  circle  target  by  an  equivalent  Gaussian  target  of  density 
exp  T  -  (X2  -r  Y*)/T*  ] 

L  J 

at  the  point  (X.  YX  the  formulas  become  simpler.  Thus  we  ftno 


? 

F 


(x.  y)  =  T  a2/ (a2  *  T2)  1  exp  [  -(x2  +  y2)/(a2  +  T2) 


=  a2/(a2  +  T2  + 


2ff2) 
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o  SF0  -  a2/ (a2  +  T2) 
Fq>  a2/ (a2  +  T2) 


SALVOS 

So  far  we  have  considered  the  probability  of  damage  with  a  single  weapon  only. 

If  several  weapons  are  dropped  separately  in  such  a  way  that  the  probabilities  are 
independent,  we  can  combine  the  single-shot  probabilities  by  well  known  methods. 

An  example  of  partial  correlation  is  that  of  a  salvo  of  weapon's,  that  is,  a  group  of 
weapons  released  simultaneously  oi  nearly  simultaneously.  In  this  case  the  weapons 
in  the  salvo  have  a  particular  distribution  which  is  ooi  the  same  in  all  respects  as 
the  distribution  that  would  be  obtained  on  another  salvo.  A  characteristic  of  the 
distribution  that  might  vary  from  salvo  to  salvo  is  the  point  about  which  the  individual 
weapons  in  the  salvo  are  grouped,  such  as  the  mean  point  of  impact  (MPI).  The 
weapons  ol  a  particular  salvo  will  have  a  particular  MPI  but  this  MPT  may  vary  from 
salvo  to  salvo.  In  focusing  attention  on  the  MPI  we  do  not  intend  te  imply  that  this 
is  the  only  characteristic  of  the  distribution  ol  weapons  iz  a  salvo  that  might  vary 
from  salvo  to  salvo,  la  the  formulation  below,  however,  we  shall  consider  only  the 
MPI. 

Let  f(x,  y  i  »i,  v)  be  the  joint  density  function  of  the  point  (x.  j)  of  (ktssztue  of 
a  weapon  in  a  salvo  that  has  its  MR  at  the  point  (u.  v).  By  the  methods  discussed 
previously  we  can  find  the  single-shot  probability,  given  chat  the  MR  is  (a,  v). 

Let  this  be  P=  P(o,  v) and  let  g(u,  v)  be  the  joint  density  hactssa  of  s  aad  t.  Lk 
r  -[?]  be  any  functional  of  ?  of  interest.  Then  the  average  raiae  of  tads  qrascetcr  is 

x 

f*  jT  -* 

\  \  F  [  ?(«,  v)  \  g(u.  v)  dr 

o'  J  t  J 

■* 

For  issuance  we  might  wars  m>  fiaf  the  procabvlsr/  of  sxrag  at  leas;  coe  daraagseg 
q«  vttli  a  salvo.  1;  the  probssbslsses  of  damage  by  weapoos  ss  a  salvo  are  asdegeaseet. 
except  for  the  coeguks  MR,  and  a  weapons  are  used  m  a  salvo,  thee  the  qrara rtsy 
of  interest  is 

f  e?3  =  i  -  a  -  p)3  =  £  t- 1)1*1  (?)  f1  <*.  *)  . 

i=l 

From  this  we  can  readily  find  the  probah&iry  sf  scornag  at  feast  ose  damagseg  fes: 
by  atrgracsg  the  product  of  this  San coos  of  cs  and  r  and  gfo.  v)  over  the  fa.  v) 
plane. 


p<e  F„> 


fi¬ 

le 


[f,>2  +  tV] 


y  o  ■» 

(a  +  i  )/2  <rA 
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We  will  treat  in  full  one  particular  example,  perhaps  the  simplest  example 
possible.  Assume  that  the  conditional  damage  function  is  circular  Gaussian  with 
equivalent  radius  a.  Assume  that  the  distribution  of  shots  in  a  salvo  is  circular 
normal  of  standard  deviation  s  and  is  centered  at  (e,  v).  Also,  assume  that  the 
distribution  of  the  MP!  is  circular  normal  with  standard  deviation  j  and  is  centered 
at  the  target.  Then  from  die  solution  to  the  problem  concerning  the  circular  normal 
distribution  of  shots  centered  at  a  point  h  units  from  the  target,  we  have 


P(u,  v)  =  A  exp 


*  <u 


2221  22  2 
+  v  )A/a  |,  A  =  a  / (a  +  2s  ) 


A  A 

T> g (u,  *1  dc  dv  =  Al/(i  *  i/BX  B  =  (a2  +  2s2V292 

~x 

Toe  latter  result  is  obtained  by  straightforward  integration-.  Finally,  the  prob¬ 
ability  of  scorisg  at  least  one  damaging  hit  cas  be  written  as 


(p(> 


1)  =  £  (*!)*  1  (f )  AVO  "  */»> 


where 


A  =  a2;**2  *  2s2) .  3  =  (a2  +  2s2y2  <r2 

The  ptobabeucv  c*  at  least  1  damaging  htt  n  this  ts  ktt'1*  »Titnf  by 
assnanaag  that  F  is  so  small  shat  we  caa  estinue  the  pnabnWlgy  sf  at  lease  1 
deicagheg  sue  when  she  MR  is  {is.  v)  aj 

F(F)~  1  -  «?  {  sPV  P=FC«.v) 

Using  this  approximation  «e  can  express  the  average  uue  n  tne  form 

* 

(?{  -  i)  -  I  -  (1/  <r  j  exp  -r2/2  s2  -  xA  exp  {-t2 A/z~)>  r  dr 


Letnag  s  =  xA  exp  (-r2A/a2)  we 

sA 


<P{2r  -  3(stA)'B  f  e^t5*1  dt  =  1  -  (xA)  8  p(»*  l)p(xA/l*3, 3-i) 

J  IB 
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where  |  (B  +  1 )  is  the  complete  gamma  function  and  (u,  p)  is  the  incomplete 
gamma  function  tabulated  in  reference  (b). 

The  above  example  indicates  the  geoeral  procedure  is  problems  of  this  kind. 

The  reader  may  wish  to  compare  the  probability  obtained  is  the  above  example 
with  that  obtained  by  first  finding  the  over-all  distribution  of  weapons,  then  deter¬ 
mining  the  probability  of  damage  with  a  single  weapon,  and  from  this  deriving  the 
probability  of  at  least  one  damaging  hit  wish  n  weapons  which  are  considered 
tndeperxleae.  The  probability  of  scoring  at  least  one  damaging  hit  by  the  latter 
method  should  be  greater  than  the  probability  found  by  the  former  (and  correct) 
method.  Weapons  fired  is  salvos  constitute  one  of  the  simplest  examples  of 
correlation. 

SEQUENCE  OF  SHOTS 

ts  asst  problems  two  or  pore  weaspaes  are  fired  oe  a  single  ran  or  e*wa  io 
soefe  a  way  that  tbe  Eadrvidtsal  probabthaes  are  not  sadepesdeet  aac  there  are  ssasy 
cfcaracteriscics  of  the  disenbessoe  that  change  from  shot  to  shoe.  Examples  are: 
5i£rface-to-air  aaaircrafe  fire  bv  guas.  rockets  aad  gaaied  pasties:  asr-t o-siz 
eegagepeses  with  gsss.  rockets,  'ssd  massifes;  stack  booh-ag  aad  sater»a*o«oeser 
boenbtag  br  aircraft.  etc.  ti  fact,  most  of  the  rertltstsc  siroaoacs  piolvisg  ai?re 
thas  ooe  weape c  are  :f  this  type.  L'soalry  the  co-rreiasso*  everts,  Noth  aaco-cor- 
reEatsre  aai  are  mo  large  to  be  Egoorsc.  Ncpe-ro'Oi  eff  MTS  a**e  been 

-rr^r-ft*  o  stride  procfatES  of  sort.  Most  of  the  s»o&EgBO<£>  recjrrre  shat  the  process 
be  sfasaocaanr.  3tet  tfw  process  is  t&suaSly  tec  statBoairy  sa  the  pcobfieaos  eccciasEeretL 

A  of  shots  that  forms  a  ssspse  Markoff  cfeace  probably  ss  the  «r=o6esc 

-Tpe.  Wsfe  a  suSnc»ea«iilj.  chars  E&rs  cethad  is  adftyace  to  treat  the  *sEo-»r* 
refiaamc  effects.  Th:  progress  of  the  Martaff  dbaaa.  mrihsdng;  the  Sseesnag  law  of 
sssasl  jussher  5  Ctthe  rmafogtae  if  the  ^season  (fejcmhiciijfia.  as  treated  ex  references 
4rJ  aad  jf).  There  seeas*  :s  be  a>  easy  sacehod  of  ireacisg  Lexiaa  effects,  except 
by  raArsg  the  wetghad  average  of  the  profiaabetoses  for  the  separate  pocarlaCDoea. 


J»ic 

i_  S7E5NKAXDT 
director 

Operaaocs  Ewlsssms  Group 


Suctsanad  by: 


J.  «-  EXMBlE 
Oaeranccs  S^abasaoc  Grorp 
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APPENDIX  A 

A  SEQUENCE  OF  CONDITIONAL- DAMAGE  FUNCTIONS 

The  conditional  damage  function  is  a  function  which  gives  the  probability  that 
a  target  which  is  located  at  a  particular  point  and  bearing  from  the  burst  point  will 
suffer  at  least  the  stated  degree  of  damage.  This  probability  is  a  function  of  dis¬ 
tance  r  of  the  center  of  the  target  from  ground  zero,  height  (depth)  of  burst, 
bearing  of  ground  zero  relative  to  the  target,  heading,  damage  category,  target  vul¬ 
nerability,  and  weapon  yield. 


It  is  evident  that  a  function  which  accurately  describes  the  effects  of  all  these 
variables  will  be  complicated.  First,  let  the  last  three  variables  --  damage  category, 
target  vulnerability,  and  weapon  yield  --  have  assigned  values.  Then  the  conditional 
damage  probability  is  a  function  of  the  three  coordinate  variables  --  horizontal 
distance  r,  height  (depth),  and  bearing.  First  average  over  the  distribution  of 
these  variables  for  a  selected  aiming  point  (sometimes  called  the  intended  ground 
zero,  IGZ)  and  intended  height  of  burst. 

The  most  important  coordinate  variable  is  horizontal  distance  r.  For  sim¬ 
plicity  the  conditional  probability  will  be  written  explicitly  as  a  function  of  this 
variable  only.  It  is  understood  that  r  particular  height,  based  on  weapon  yield, 
effect  desired,  dud  probabilities,  etc.,  has  been  set  in  the  fuze.  It  is  assumed 
that  errors  in  height  produce  a  negligible  effect  on  the  probability  of  damage,  or 
that  an  average  over  the  corresponding  error  distribution  wiii  be  taken.  Such  an 
average  can  be  taken  at  this  point  if  the  height  error  is  independent  of  errors  in  the 
horizontal  plane,  which  is  true  for  most  delivery  and  fuzing  systems  in  use.  How¬ 
ever,  for  air-burst  fuzes  it  is  found  to  be  convenient  to  defer  this  average  until 
later. 

It  is  assumed  that  changes  in  bearing  produce  negligible  changes  in  damage 
probability.  For  some  targets,  notably  ships  and  aircraft,  this  is  not  true.  How¬ 
ever,  in  taking  the  average  over  the  distribution  of  weapons,  these  changes  with 
bearing  can  usually  be  ignored  with  negligible  error. 


When  the  conditional  damage  probability  is  written  as  a  function  of  r,  it  will 
be  assumed  that  the  effects  of  changes  in  damage  category,  target  vulnerability, 
and  weapon  yield  can  be  described  adequately  by  changes  in  paxameters  in  this 
function,  without  changing  the  form  of  the  function. 


Let  p  (r)  be  the  average  probability  of  damage  at  horizontal  distance  r  for  a 
given  damage  category,  target  type,  and  weapon  yield.  Several  functions  which  have 
been  used  are  tfce  following: 


(a)  Definite- range  or  cookie-cutter  function: 


P0<r>  = 


1 


if  r  s  R 

if  r  >  R 


(A-l) 


where  R,  the  damage  radius,  is  a  parameter  which  depends  upon 
damage  category,  target  type,  and  weapon  yield. 


t 


A-l 


(b)  Gaussian  function; 

Pj(r)  =  exp  (-  r2^2) 

where  "a^  “  also  is  a  parameter  which  depends  o*  tie  damage 

variables.  A  facto;  of  2  is  sometimes  used  in  the  dewosciastnr  of  tbe 
exponent  to  simplify  later  formulas;  this  should  be  kept  is  mind 
when  comparing  values  of  parameters. 

The  definite -range  function  appears  to  be  reasonable  if  ail  the  r&riaiaes  except 
r  have  f  jced  values.  This  requires  that  the  damage  category  be  defined  precisely, 
a  particular  target  of  known  vul  >er  ability  be  considered,  the  weapon  yield  be  known 
exactly,  and  die  height  and  bearing  of  the  burst  position  be  fixed.  As  cxpiais^d 
above,  p(r)  should  be  an  average  over  some  of  these  variables  --  particular!? 
variations  in  damage  which  would  be  included  in  a  given  damage  category,  varsctxaes 
in  target  vulnerability  among  targets  of  a  given  type,  aad  variances  ca  atigja.  c£ 
burst  unless  the  average  over  the  height  distribution  is  to  be  takes  later. 

For  any  conditional  damage  function  p(r  j.  let  g(R)  be  the  carrtspotthFj  dessty 

function  of  tbe  damage  radius  R.  Write  p  (r)  ta  the  form  p  (r,  R)  to  cksplry  the 

o  *  o 

role  of  R.  Then 


P<r)  =  J 
o 


Po(r.R)  g  (R)  dR 


ar 

\  g(R)dR 

V 

r 


(A* 


•n 


is  the  distribution  function  of  the  damage  radios.  Hence,  the  probfcaa  of  finding 
the  probability  function  p(r)  is  equivalent  to  that  of  finding  the  (hstrihatioc  fsactisc 
of  the  damage  radius  R. 

For  p(r)  equal  to  the  Gaussian  function  p  (r)  the  correscowding  density  fsoctsas* 
of  R  is  1 


2R  2  2 

gjCR)  =  -  pL'  (R)  =  — y  exp  (-  R  /a') 

al 

An  objection  to  ibis  ‘density  function  is  that  it  gives  too  much  weight  to  very 
small  and  very  large  values  of  R.  Asjkill  be  shown  later,  only  65  percent  of  the 
distribution  lies  between  0. 5lC  audi.5R,  where  R  is  the  average  damage  radius. 

To  get  a  higher  concentration  about  the  average  we  can  multiply  the  expooessial 
factor  by  a  power  of  R  greater  than  oue.  For  example,  if 

3 

g2(R)  =  exP  (~  2R2/a2)  , 

£2 


A-2 


i-l 


JL*  X  2 


*  «5T3«aoI 

«fc 

E  'tt  arfartf  Jn 
s 


1  I-  i 

Fir  dr  Cfctacsurr  iacag  dr 

E  5  *»  1  a**  1  ;  t_C  x. 
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FIG.  A-l:  DENSITY  FUNCTIONS  OF  DAMAGB  RADIUS  R 


PROBABILITY  OF  DAMAGE  AS  FUNCTION  OF  HORIZONTAL  RANGE 


2  _ 2 

Heace,  the  variance  from  tilt  mean  is  R  -  R  .  The  standard  deviation  from 
_  _  e 

the  mean  is  0. 52R ,  0. 36R,  and  0. 29R  for  (r),  p^  (r),  and  (r)  respectively. 

If  we  know  R  and  the  standard  deviation  of  R  from  R  ,  the  appropriate 
function  can  be  chosen  from  the  sequence  of  functions.  If  only  R  is  known,  how 
sensitive  is  the  probability  of  damage  to  the  choice  of  conditional  damage  function? 

Let  (x,  y)  be  the  rectangular  coordinates  of  ground  zero  relative  tc  the  center 
of  the  target.  Assume  that  x  and  y  are  independent  and  normally  distributed  with 
the  same  standard  deviation  o  and  means  h  and  h  respectively.  (Unequal 

^  y 

standard  deviations  make  the  formulas  more  complicated  but  do  not  increase  the 
difficulties  of  integration. )  The  density  function  of  x  and  y  is 


f  (x,  y)  = 


["<x‘ 

\a  I 


hx)  +  (y  -  h 


2*  O' 


Then  the  average  probability  of  damage  for  the  n—  function  is 


’°=iT Pn 


(x,  y)  f  (x,  y)  dxdy 


where  p^(x,  y)  is  the  function  given  in  equation  (A- 3)  when  r  is  replaced  by 

t  2  2. 

(X  +  y  ) . 

2  2 

Expanding  the  powers  of  (x  +  y  )  and  completing  the  squares  in  the  exponent, 
it  is  evident  that  the  integral  in  Pq  can  be  computed  for  any  n.  We  omit  details 

and  give  the  result  after  the  integration  is  completed.  At  this  stage  we  have 


-H 

p  = e  a  -  aJ  S, 


n'  n 


(A-4) 


where 


n-1  n-1 
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i  =  L  La  — p -  S  P 

k=o  j=k  2J  •  j  :  Jk  n 


(A- 5) 
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(I) 


(A- 6) 


A-7 


He  factor  5^  nr  Jani  %  to 

sjk  *  if  <i>  *• 

******  das  uco  (A- 5)  we  hm  dat 


»-i 

z 

k=o 

£ 

j=* 

(IK 

V 

HET 

«-l 

fr:): 

..  k 

I 

k=o 

*£.  U»> 

t  S 

(V 

i*cre  («#)  **  tie  k  -  dttivaime  td 


*x  =e'Hl(l  -  «  p 

.U 

p2 =  e  (A  - a  2)  [  a  +  «  2)  +  a  -  o  j)  h2  ] 


A-8 


I 


TV  1 3C-'  yae  ttm. 
z  -  «— c*nuar-raay 

**c  1  r  it  uew  * 

dae  s€  F  .  dr 


kt.i  -  d*«  «-y  ter  *rt  dcuari  wtt  x  =  :  jOwcag)  ** 
TV  cairt  ;f*  F  »_  m  F  .-aaetf  ter  —at  «ri—  «f 

.  M 

F’C  «suT  Kxian  aft  r  At  ni*  s£  ?.  a  f^rwm ~  d«a 
Meat  .»  rrat  3:r  nuua  W.  e  .  TV  vtut  ai\  t 

s  ».!ic  ;  r  i  ..I  fa:  k  #  :  l.  *a£  dMt  If  fa: 

«res  F  as2  ?  s  »  m  *s  >r  yrggg—  tt  tacratat  run 


I  m£ 

S  iff  f 


siaasaa 


FIG.  A-3:  CONSTANT  PROBABILITY  CURVES  FOR  ?2 
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